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We discuss preliminary results from our quenched light hadron mass calculation on a 48^ x 64 lattice at the 
coupling of /3 = 6.5. Staggered quarks with masses of niq = 0.01,0.005,0.0025 and 0.00125 are used. 



1. MOTIVATION 

We are calculating light hadron masses 
in quenched lattice quantum chromodynamics 
(QCD) with staggered quarks. In a series of works 
[1-3], one of the authors (SK) was involved in a 
study of the systematic errors arising from the fi- 
nite lattice spacing, finite volume, and finite quark 
mass. An important suggestion from those works 
is that a combination of the spacelike lattice size 
of L = 2.4 fm and the coupling of /3 = 6.5 should 
be good enough for extrapolations to remove the 
errors from finite volume and finite lattice spa- 
cing. The errors from finite quark mass is more 
subtle. Sharpe [4] and Bernard and Golterman 
[5] suggest that the chiral limit of quenched QCD 
differs from that of full QCD. In ref. [3], SK and 
Sinclair argued that indeed the chiral behavior 
of quenched QCD is different from the full the- 
ory in regard of the suggested "quenched chiral 
log" term. In those works [1-3], however, the 
three sources of errors were studied separately 
due to the hardware limitation. In contrast the 
present work simultaneously addresses the finite- 
volume and finite-lattice-spacing effects. It uses 
a VPP500/30 vector-parallel supercomputer at 
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RIKEN [6]. We plan to investigate the chiral limit 
of quenched QCD in more detail, particularly in 
regard of a suggestion that the finite volume ef- 
fect compounds the chiral limit [7] . In the current 
study the chiral limit should be approached more 
reliably than in the stronger coupling cases since 
the fiavor symmetry violation is small. 

2. SIMULATION CHARACTERISTICS 

The 48^ X 64 lattice volume and the coupling of 
/3 = 6.5 are selected for the reasons discussed in 
the above. The spacelike lattice size of 48 makes a 
24-node partition of the 30-node VPP500/30 con- 
venient. As was already reported [6], each node of 
this computer has a peak speed of 1.6 GFLOPS, 
256 Mbytes of SRAM memory and an ability to 
communicate with any of the other nodes at 800 
Mbytes/s peak. Through coding in a modified 
fortran we obtain more than 1 GFLOPS sustained 
per node for the current work. 

In generating the gauge configuration we use a 
Metropolis update followed by an over-relaxation, 
as in the earlier works [1-3]. The separation 
between hadron-mass calculations is 1000 such 
pairs of Metropolis and over-relaxation sweeps, 
and take about 3 hours in total including the 
necessary disk access. Since decorrelation took 
about 500 such sweeps for the 32^ x 64 lattice 
at /3 = 6.5, the 1000-sweep separation should be 
good enough. We have not seen any sign that 
the contrary is true. All the configurations used 
for the hadron-mass calculations, almost 2 Gbytes 
each, are stored in a tape archive. This will en- 
able us to study hadrons with strangeness, charm 
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and bottom in the near future. 

Before calculating the propagators, we calculate 
chiral condensate {i>il>) for each masses using bi- 
linear noise scheme with single noise vector. 

The gauge field is fixed to Coulomb gauge by an 
over-relaxation method prior to quark-propagator 
calculation. For estimating quark propagators we 
use the standard conjugate gradient (CG) method 
using fixed 48^-volume corner walls and even 
point walls as sources and a point sink. Staggered 
quarks with masses of = 0.01,0.005,0.0025, 
and 0.00125 are used. The physical wall size 
should be about (2.4fm)^ and the physical quark 
mass should range from about 5 MeV to 40 MeV 
since the lattice cutoff is ~ 4 GeV at j3 = 6.5 
according to the scale set by the p meson mass [1]. 
The convergence criterion for the CG iterations is 
the norm of the residual vector being equal to or 
smaller than 1.0 x 10~^. It takes from about 1000 
(for niq = 0.01) to 8000 (for rUg = 0.00125) CG 
iterations till convergence. The calculations from 
the chiral condensate and gauge fixing to hadron 
propagators take about 6 hours in total for each 
gauge configuration. All the quark propagators, 
almost 1 Gbytes each, are also stored in the tape 
archive. 

After the hadron propagators are calculated, we 
fit them using the same procedure as in the previ- 
ous studies [1-3]: CERN library MINUIT is used 
for minimization of the correlated function and 
the error bar quoted is for one unit change in ■ 
At the time of the Melbourne meeting we repor- 
ted the results from 30 gauge configurations. Here 
we choose to report the results from 50 configur- 
ations, an increased statistics in the month fol- 
lowing the meeting, and using the full covariance 
matrix. 

3. RESULT AND DISCUSSION 

In Table 1 and 2 we summarize, respectively, 
the chiral condensate and light hadron masses 
for each of the four quark mass values calcu- 
lated from the 50 gauge configurations. Note that 
all the results presented here are still prelimin- 
ary. Complete analysis should wait till we ac- 
cumulate twice more statistics or 100 configura- 
tions. In the tables we also compare the current 



Table 1 



Chiral condensate {tptp) at /3 = 6.5. 



niqa 


48^ X 64 


32^ X 64 


0.01 


0.04467(3) 


0.04466(4) 


0.005 


0.02382(4) 


0.02387(7) 


0.0025 


0.01337(6) 


0.01348(9) 


0.00125 


0.00808(8) 





50-configuration results on the 48^ x 64 volume 
with those of the 100-configuration ones on the 
32^ X 64 lattice [1]. From this comparison it is 
clear that the current preliminary results from the 
48^ volume are consistent with the established res- 
ults from the 32^ one. On the other hand we still 
do not understand why the current results tend to 
give smaller error bars despite their lower statist- 
ics: perhaps through smaller fiuctuations due to 
the larger lattice volume, or perhaps the x^ fit- 
ting procedure is under-estimating the error. A 
different error estimate such as those based on 
the Jack-knife method would be helpful, and is 
planned after we accumulate more data. 

It is also worthy to note that the current pion 
effective mass, plotted in Figure 1, gives better 
plateau than that from the previous 32^-volume 
calculations though the noise sets in early in time 
coordinate due to lower statistics. This suggests 
that the (2.4fm)^ wall is indeed optimal. Interest- 
ingly, in the previous calculations [3] at /3 = 6.0, 
effective mass plot for 24^ x 64 also shows better 
plateau than that for 32^ x 64 or 16^ x 64. 

Figure 2 gives the Edinburgh plot. The errors 
in the figure are obtained by assuming that the rel- 
ative error in each quantity is independent of each 
other. In the plot we used the nucleon mass from 
the "even-point-wall" source because the quality 
of fitting is better than that from the "corner- wall" 
source. The other quantities are from the "corner- 
wall" source. In contrast to the same plot for 
the 32^ volume, all the points, except for the left- 
most point for the lightest quark mass, lie below 
the guiding lines. This suggests that ground-state 
nucleon is less squeezed on the 48^ volume than 
on the 32^. 

At the lightest quark mass of = 0.00125, we 
now see the ratio m-^ /mp as low as ~ 0.3. Since 
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Table 2 



TTarlmn inaQQpQ 


at /? — fi 5 on 48^ > 


\ VJt: 1 Cli Lj Lj 1 <3 . 






TJlq a 


particle 


fitting TRllgG 111 t 


ihrss 


mass trom 32 




TT 


4-17 


0.1630(6) 


0.1592(10) 






6-25 


0.1605(7) 


0.1611(12) 




P 


8-17 


0.2489(15) 


0.2417(24) 




P2 


7-17 


0.2483(13) 


0.2416(36) 


0.01 


ai 


7-17 


0.3489(30) 


0.3481(49) 




hi 


10-17 


0.3406(80) 


0.3498(81) 




a 


6-20 


0.3062(42) 


0.3197(62) 




Ni 


10-18 


0.3678(25) 


0.3761(34) 






3-18 


0.3444(17) 


0.3522(29) 




A 


7-16 


0.4141(26) 


0.4193(48) 




TT 


6-18 


0.1179(10) 


0.1133(13) 






6-20 


0.1172(12) 


0.1162(18) 




P 


8-20 


0.2285(23) 


0.2179(42) 




P2 


7-17 


0.2312(23) 


0.2208(59) 


0.005 


ai 


5-17 


0.3261(33) 


0.3296(80) 




hi 


8-20 


0.3251(73) 


0.3348(186) 




a 


6-20 


0.3044(94) 


0.3455(210) 




Ni 


10-20 


0.3187(44) 


0.3416(58) 






12-31 


0.3013(50) 


0.3106(54) 




A 


5-29 


0.3989(29) 


0.3934(63) 




TT 


6-18 


0.0855(12) 


0.0821(17) 






6-20 


0.0884(19) 


0.0893(34) 




P 


8-20 


0.2186(37) 


0.2043(65) 




P2 


4-18 


0.2317(26) 


0.2049(109) 


0.0025 


ai 


6-28 


0.3076(38) 


0.3063(67) 




hi 


5-18 


0.3228(72) 


0.3323(468) 




a 


6-20 


0.3233(200) 






Ni 


8-19 


0.3045(80) 


0.3287(115) 




#2 


11-29 


0.2614(93) 


0.2859(88) 




A 


4-22 


0.4007(42) 


0.3683(157) 




TT 


6-17 


0.0633(15) 






T^2 


4-13 


0.0608(38) 






P 


4-18 


0.2214(39) 






P2 


4-18 


0.2338(52) 




0.00125 


ai 


4-18 


0.3028(51) 






hi 


5-18 


0.3083(134) 






a 


3-13 


0.2735(155) 






Ni 


4-20 


0.3318(82) 






#2 


3-23 


0.2891(58) 






A 


4-22 


0.4029(52) 
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Figure 1. Pion effective mass at [3 = 6.5 on 
48^ X 64 iattice for quark masses niq = 0.01, 0.005, 
0.0025 and 0.00125. Corner-wall source. 



Figure 2. Edinburgh plot at (3 = 6.5 on 48^ x 64 
lattice for quark masses niq = 0.01, 0.005, 0.0025 
and 0.00125. The two circles are the experi- 
mental point (left) and naive quark model pre- 
diction (right). 



we have m-^L ~ 2.9, we can probably conclude 
that our lattice is big enough even down to such 
a light quark mass. Thus pushing the quark mass 
niq even smaller and reaching the region where 
m-^lmp ~ 0.2 with m-^L ~ 2.0 and m-^ ~ 0.04 
seem viable. In other words extrapolation to small 
niq in the quenched spectrum calculations may 
soon become unnecessary. 

Also the current result shows approximate fla- 
vor symmetry restoration: the masses of tt and 7r2 
agree with each other for all the calculated quark 
masses and the masses of p and p2 agree for the 
two heavier quark masses and are close to each 
other for the lighter two. 

All these attractive features of the current res- 
ults suggest that the chiral limit of our data will be 
more interesting. However, we think that the in- 
vestigation on the quenched chiral log in the pion 
mass m,r and chiral condensate {i>il>) should wait 
for a still higher statistics of at least 100 gauge 
configurations. We do not think the current stat- 
istics is enough for that purpose because a) each 
point in the pion effective mass plot shows large 
fiuctuation despite the fact that the error estim- 



ated by the full covariance matrix is small, and 
b) the /9-meson mass data also imply large fiuc- 
tuation as nip at niq = 0.00125 is similar to, and 
not smaller than, that at niq = 0.0025. 
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